Exotic hadron production in quark combination model 
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The philosophy on production of exotic hadrons (multi-quark states) in the framework of quark 
combination model is investigated, taking /o(980) as an example. The production rate and pr 
spectra of /o (980) considered as (ss) or (sqsq) respectively, are calculated and compared in Au+Au 
collisions at y^ijvjv ~ 2006*61^. The unitarity of various combination models, when open for exotic 
hadron production, is addressed. 
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I. INTRODUCTION 



The basic quanta of Quantum ChromoDynamics 
(QCD), quarks and gluons, are confined in their bound 
states, hadrons. In any high energy interaction, the pro- 
duced colour-singlet (CS) (anti)quark system eventually 
transits to various hadron states with the total probabil- 
ity exactly 1: 



h 



< h\U\q > P =< q\U+U\q >= 1. 



(1) 



Here we introduce the unitary time-evolution operator 
U to describe the hadronization process. For the quark 
state \q > and the corresponding hadron state \h >, the 
matrix element Uhq —< h\U\q > describes the transition 
amplitude. Uhq is determined by QCD but beyond the 
present approach of calculation. This leaves the space 
for various hadronization models to mimic this transition 
process. As a matter of fact from experiments. 



E I 

h=B,B,M 



< h\U\q > 



1 - £, e -> 0^ 



(2) 



here B,B,M denote baryon, antibaryon and meson re- 
spectively. 

Naively from the group theory, colour confinement 
seems not so restrict as Eq. The CS state, 

i.e., the invariant, totally antisymmetric representa- 
tion of the S'[/c(3) group, requires at least one quark 
and one antiquark, or three (anti)quarks (These are 
just the constituent/ valence quark numbers for mesons 
and baryons, respectively), but more (anti)quarks can 
also construct this representation, hence possibly to 
form a CS "hadron" . To name some possibilities, two 
quark- antiquark pairs, a quark-antiquark pair with 3 
(anti)quarks, six (anti)quarks, etc., are to be called exotic 
hadrons in this paper 30]. Until now, no experiment can 
definitely show e in Eq. ^ is exactly or a small but 
non-vanishing number. If definitely e = 0, there must 
be underlying properties of QCD which we still not very 
familiar. Even s is not vanishing, its smallness, definitely 
confirmed by experiments and shown in Eq. ([2]), also pro- 
vides interesting challenges, especially on hadronization 



models. The small production rate of a special kind of 
exotic hadron seems easy to be adopted. However, tak- 
ing into account so many possibilities to construct the 
CS representations by various numbers of (anti)quarks, 
that the total sum of them is still quite small, is very 
nontrivial as property of QCD and even nontrivial for a 
hadronization model to reproduce. 

More concretely, we investigate the quark combination 
models when open for the production of multi-quark 
states. At first sight, the quark combination model is the 
most feasible in calculation, to just allow the desirable 
number of (anti)quarks to combine together. But to get 
the total production probability of all the exotic hadrons 
fulfilling Eq. ([2]) with a universal combination rule is 
nontrivial. In this paper, we investigate the production of 
exotic hadrons via a quark combination model proposed 
by the Shandong Group (SDQCM) 0, i, [s, 6, 1, 8]. We 
demonstrate that to treat the production process of all 
multi-quark states as hadron molecule formation could 
be a practical and self-consistent way, respecting Eqs. 
([11 [2|) , the unitarity of hadronization models (see also |9| , 
and refs. therein). 

Of all the "on market" combination models, SDQCM is 
unique for its combination rule wisely designed, so that 
mesons and baryons exhaust the probability of all the 
fates of the (anti)quarks in a CS system. It has been real- 
ized in Monte-Carlo programmes and tested against data 
from experiments in high energy e~^e~ annihilation and 
pp collisions, and recently it has been successfully extrap- 
olated to ultra- relativistic heavy ion collisions fiol. [ill. . 
refiecting the universality of the hadronization mecha- 
nism. Most recently, application of the combination rule 
to the open charm and bottom sector for Relativistic 
Heavy Ion Collider (RHIC) experiments, without any 
more parameters for the hadronization, further demon- 
strates its validity and provides opportunities against 
more critical tests [Tsj . 

The main idea is to line up Nq (anti)quarks in a one- 
dimensional order in phase space, e.g., in rapidity, and 
then let them combine into initial hadrons one by one 
following a combination rule [1, 0, H, B S Hi) three 
(anti)quarks or a quark-antiquark pair in the neighbor- 
hood can form a (anti)baryon or a meson, respectively. 
The cases of more numbers (> 3) of (anti)quarks to com- 
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bine into hadrons automatically disappear. The inclu- 
sive cross section is proportional to the product of the 
quark number densities, leading to B/M enhancement 
in a region (which is in fact one of the motivations for 
the burst of combination models at RHIC). Without this 
special combination rule, more-quark state could be even 
enhanced [§] . This will conclude a nonsense result that a 
system of (anti)quarks could combine into a giant "quark 
ball" , with large number of constituent quarks. So, one 
of the key properties of combination models is whether 
a model can shut off the possibility that exotic number 
of (anti)quarks combined into cluster with a large prob- 
ability. 

To introduce the small probability e in Eq. ([2]) to 
SDQCM, similar as the framework to calculate the res- 
onance production of heavy quark bound states p^ . 
we consider the exotic hadron production as hadron 
molecule production and project the free mesons and/or 
baryons (anti-baryons) states onto their bound states in 
terms of the wave functions. In this paper we first clarify 
that it does not mean that, in the bound state, quarks 
from each hadron should keep in CS respectively, but the 
colour interactions can transit the whole molecule into a 
"real" exotic hadron by some probability, which will be 
introduced in the following. This is especially reasonable 
for the circumstance like RHIC, where the large num- 
ber of hadrons and large bulk of thermalized area allow 
interactions between hadrons lasting for a long time. 

In the following we first discuss the colour state and 
"definition" of quark number of the exotic hadrons (Sec. 
II). Then, taking /o(980) as a working example, we de- 
scribe the calculations of the production of the exotic 
hadrons within the quark combination model and discuss 
the results (Sec. III). Sec. IV is the conclusion. 

II. COLOUR STATE AND QUARK NUMBER 
OF THE EXOTIC HADRONS 

A. Colour state in an exotic hadron 

All kinds of Exotic hadrons have one common prop- 
erty, which is that the (anti)quarks can be grouped 
into several clusters, with each cluster possibly in CS. 
Hence, could the exotic hadrons just be meson and/or 
baryon molecules? However, the ways of grouping these 
(anti)quarks are not unique, as it is simply known from 
group theory that the reduction ways for a direct product 
of several representations are not unique. Furthermore, 
these clusters need not necessarily be in CS respectively, 
since only requirement is the whole set of these clusters in 
CS. For example, the system 91(729394 (the constituents 
of a "tetraquark" ) can be decomposed/clustered in the 
following ways: 

(9193)3 ® (9"294)3 -> 1 (3) 
(9192)1 or 8 ® (9394)1 or 8 1 (4) 



Here we just mention that such group theory analysis is 
applicable to the quark states as well as the quark field 
operators [l^. In the above example, only the second 
case, when these two 99 pairs are in CS respectively, it 
seems possible to be considered as a hadron molecule. 
But dynamically, the colour interactions in the system 
via exchanging gluons can change the colour state of each 
separate cluster, so each kind of grouping/reduction way 
seems no special physical meaning. Such an ambiguity, 
which has been considered in many hadronization and de- 
cay processes as "colour recombination/rearrangement" 
[3 Ell , obstacles the possibility to consider the ex- 
otic hadron in a unique and uniform way, while leads to 
the possibility of introducing some phenomenological du- 
ality. Namely, even we consider the production of exotic 
hadron as "hadron molecule" formation, the subsequent 
colour interactions in the system can eventually transit 
this "molecule" into a "real" exotic hadron, at least by 
some probability. 



B. How to count the quark number in an exotic 
hadron 

Some kinds of exotic hadrons have exotic quantum 
numbcr(s), e.g., one kind pentaquark (qqqsq), has +1 
baryon number but -fl strange number. If a hadron 
with such quantum numbers is experimentally confirmed, 
one seems to have to introduce five valence (anti)quarks. 
However, in many other cases, there exist parallel ex- 
planations because of non-exotic quantum number of a 
certain hadron. The tetraquark or four-quark state dis- 
cussed in this paper is an example. One of the candidates 
is /o(980). It is considered as an orbit-excited I = I reg- 
ular meson (ss) [l9j , but argued possibly to be a four- 
quark state {sqsq) by others |20j . So one naturally raises 
the question relating with this ambiguity, how on earth 
can one count the number of quarks in a hadron? Even, 
what is the meaning of "number of quarks in a hadron" ? 
This may be one of the most ambiguities in the physical 
picture of the exotic hadrons, because of lack of complete 
understanding of the confinement property of QCD. It is 
not clear what the quantum field theory definition of the 
"constituent quark" is. As a consequence, it is ambigu- 
ous how to "count" the quark number in a hadron. Here 
we just state the different pictures of a proton, one is par- 
ton model bursting from the deeply inelastic scattering 
and other high energy interaction processes, the other is 
the "static" quark model, corresponding to the properties 
of a proton at rest. There is no satisfactory, especially 
quantitatively relation between these two pictures. Even 
the consideration of higher Fock states or pair excitations 
[2l| can not remove the gap. 

It is well-known that factorization theorem confirms 
the parton fragmentation picture to describe the hadron 
production [22] . But this partonic picture is only valid in 
inclusive processes with hard interactions involved. For 
the low pt particles, e.g, those from "fragmentation" of 
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the hadron remnant in hadron-hadron interactions, or in 
a very complex circumstance like in heavy ion collisions, 
this picture faces both challenges from theoretical as well 
as experimental aspects. It is difficult to prove the fac- 
torization theorem for these complex cases. Furthermore, 
the RHIC data expose some properties difficult to be un- 
derstood from fragmentation picture, such as the high 



p/tt ratio at intermediate transverse momenta 23 1 and 

21, but 



the quark number scaling of hadron elliptic flows 
these can be explained by coalescence or (re) combination 
models [111, [2^ . Combination models include the follow- 
ing picture: 1) the production of quarks, which is con- 
sidered as the "dressed" quarks, i.e., constituent/ valence 
quarks; 2) these quarks combined to certain hadrons, 
e.g., qq combined to meson. Within these models, the 
quark number in a hadron has a model-dependent but 
clear "definition", i.e., "the number of (anti)quarks" of 
a certain hadron is the number of (anti)quarks involved 
in combining into the hadron in the production process. 
So, as suggested by [2^, one can "count" this number 
experimentally by measuring the V2 in heavy-ion colli- 
sions at RHIC, since V2 is proportional to the number of 
quarks combined to the certain hadron, and the follow- 
ing inner interactions between (anti)quarks in the hadron 
will not change this property of global movement. This 
makes the quark number an "observable" , by extrapolat- 
ing the combination picture to open for all kinds of exotic 
hadrons. Treating exotic hadron production as hadron 
molecule formation will not change this fact that V2 is 
still proportional to the total quark number in the exotic 
hadron, which is a straightforward result of associative 
law of addition. 



III. THE PRODUCTION OF /o(980) AT RHIC 

Now we apply the above discussions into an example, 
/o(980). This particle is considered as a tetraquark with 
the flavour content sqsq The traditional consider- 

ation for it as / = 1 state of ss, has been included in 
SDQCM with all the other I — 1 hadrons in the most 
recent version [28] . We also calculate that case for com- 
parison. 

As we mentioned above, the conventional SDQCM can 
calculate the free meson and baryon distributions, e.g., 
two meson distributions ^ipf^ip^ ■ Projecting onto the 
meson molecule state, we can get the bound state distri- 
bution: 



Eda 



N 



d^p 



r d^pi d^p2 

^ J El E2 

N1N2 

EiE2da^'^^^ 
d^pid^P2 



< pi,p2,Ni,N2\p,N > 



(5) 



In the above equation, the N,Ni,N2 refer to discrete 
quantum numbers, corresponding to the meson molecule 
state and the two mesons, respectively. The conservation 
relations like S^{p — pi — P2) and 5mi®N2.n are indicated 



in the projection of the state vectors. In order to get this 
factorized formulae, some interference terms are missed. 
Since the free hadron distribution is calculated by the 
Monte-Carlo programmes, which can not give the ampli- 
tude but only its square. And only this factorized form 
can be used in the Monte-Carlo programmes. 

The projection of the discrete quantum numbers such 
as flavor, isospin, angular momentum, charge and space 
parity etc., are easy to be calculated. In the case of an- 
gular momentum projecting, spin counting is assumed. 
Though the system is not assumed to be non-relativistic, 
only the lowest possible orbit angular momentum is con- 
sidered (i.e.. Since in generally, the Parity requires odd or 
even / value, in the case of even we only consider I = 0, 
while for odd, only / = f ). The detailed analysis on 
the projection of discrete quantum numbers is outlined 
in Table U 

The phase space wave function of the exotic hadron in 
terms of mesons and/or baryons is not definitely avail- 
able. One of the reasons is that, as argued above, the 
subsequent colour interactions in the "molecule" ruin a 
unique structure to be described by definite wave func- 
tion. To mimic the combination process, we use the same 
physical picture as quarks combining into regular hadrons 
in our model, i.e., near rapidity correlation [3]. This says 
that only two mesons in the neighborhood on the ra- 
pidity axis have the chance to combine together. Con- 
trary to the case of quarks, in which confinement prop- 
erty requires all quarks must be combined into certain 
hadrons, the hadrons need not necessarily to be com- 
bined to some molecules. Rather, they may have large 
probability to be free. It is natural to introduce a param- 
eter X, smaller than 1, to parameterize the probability for 
two mesons neighbored to combine into a bound state. 
So, X not only reflects the information of some specific 
exotic hadrons, but also reflects the interactions among 
the meson and (anti)baryon system. This assures the to- 
tal results respecting the fact that all the quark states 
and all the hadron states (mesons, (anti)baryons and ex- 
otic hadrons) respectively form two complete sets of bases 
for the same Hilbert space. This unitary transformation 
assures the unitarity of the quark combination model. 

Unfortunately, the unknown wave functions of all kinds 
of exotic hadrons mean x is almost a free parameter, only 
constrained by the data of mesons and baryons, as well as 
how many kinds of exotic hadrons to be considered. It is 
clear that to an extreme if we have infinite kinds of exotic 
hadrons, x should be vanishing, expecting infinite num- 
ber of vanishing variables (production rates correspond- 
ing to each certain exotic hadron) summing up to get a 
finite small result (the total production rate of all exotic 
hadrons). For demonstration, here we only consider the 
existence of four-quark states as exotic hadrons, but still 
with two choices: One case is that only /o(980)(sqs(7) ex- 
ists in the world, the other is not only /o (980) (sgsq) but 
also any other tetraquark states (isospin muti-states) to 
be allowed to produce by our model. In the latter case 
more mesons are "used up" to get the molecule, so the x 
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meson pair 


1 < Ji, J2IJ > 1^ 


1 < Ci,C2|C > {-^ 


1 < /l,/2|/ > \' 




1/9 


1 


1 


'h 'I 


1 


1 


4/9 


v,v' 


1 


1 


5/9 


V ,v 


1 


1 


4/9 


K+,K- 


1 


1 


1/2 




1 


1 


1/2 



TABLE I: The projection of discrete quantum numbers (I, 
C, P) of Ml, M2 to /o (980) (sgsg). The cases with much 
smaller values of | < pi,P2, Ni, N2\p, N > l'^ are included in 
the programme but ignored in the above list. 





X 


/„,MS(,H.,,^,,, 


/o(980)(ss) 


(a) 


0.60 


1.63 


0.68 


(6) 


0.24 


0.65 



TABLE IL The production rate of /o(980) at midrapidity 
(within one rapidity unit) for central Au+Au collisions at 
y^SNN ~ 200GeV. (a) corresponds the case only one exotic 
hadron, i.e., /o (980) (sgsg); (fe) denotes the case that we con- 
sider all the isospin multi-states. 



must be smaller than the former case. We first tune the 
Monte-Carlo programme to produce the multiplicities of 
mesons and baryons at the central value of the experi- 
mental data. The baryons are not used to produce the 
exotic hadrons, so is definitely determined by experimen- 
tal data. Then we tune x so that all the mesons used up 
should not exceed the error bar of the experimental data 
(within 5 per cent [3l|). The results are in Table [Til and 
Fig. 1. From the above discussion, it is more fair to take 
them as uplimits of the production rates. 

Recently, argued in Ref. [27|, from the production 
rate, one can extract the information of the constituent 
quark number. The authors employ the combination 
model proposed by [131, predict and compare the pro- 
duction rate of /o(980) for two cases: ss or sqsq. How- 
ever, whether this idea is practical relies on whether it is 
"model-independent" . From our calculations list in table 
im one can see different results from [23] ■ If we take into 
account the case that many other kinds of exotic hadrons 
could also exist in nature and formed from combination 
of hadrons, the production rate of /o (980) (sgsg) should 
be even smaller. So, this implies that the production 
rates depend on the mechanism of production, and are 
not straightforwardly possible to be related with the con- 
stituent quark number in the framework of various quark 
combination models. 

For the transverse momentum spectrum (Fig.l), it is 
clear the four quark state is harder, common property of 
all combination models. This is the same reason for B/M 
enhancement in mid-pT region. It is quite interesting to 
point out that, in principle the pt spectrum of such kind 
of exotic hadrons formed from combination of mesons 
and/or (anti)baryons can be fixed in our model, indepen- 
dent of the quark distributions. The free hadron spectra 
in Eq.(5) can be fixed by experimental data, which is 




FIG. 1: Hadron spectra of /o(980) at midrapidity as a func- 
tion of transverse momentum pr for central Au+Au collisions 
at ^ysNN ~ 200GeV . /o(980)(ss) (solid line) compared with 
fo{980){sqsq) (Dashed line), (a) and (b) correspond to a and 
b in table [IT] respectively. 



completely model- independent. Furthermore, x can only 
change the absolute value but not the shape of spectrum. 
So, if some exotic hadron (maybe not /o(980)) is pro- 
duced by combination of regular hadrons at RHIC, we 
can predict the shape of its spectrum without any ambi- 
guity. 



IV. CONCLUSION 

In this paper, we propose a way to calculate the 
production rates of exotic hadrons (multi-quark states) 
within the framework of SDQCM. We point out that this 
special combination model employs a specific combina- 
tion rule to shut off the possibility that more than three 
(anti) quarks can combine into some hadron. The unitar- 
ity is automatically kept but may be too restrict if there 
could be exotic hadrons, since the to-date experiments 
only assure the rareness of exotic hadron production, as 
shown in Eq. ([2]). By an analysis on the complexity of 
the colour structure of the exotic hadrons, we suggest 
that one can introduce the small section of exotic hadron 
production by treating (only) its production process as 
hadron molecule (bound state) formation. Taking the 
possible four-quark structure of /o(980) as an example, 
we discuss two special cases, i.e., only one exotic hadron 
fo{980){sqsq) exists in nature, or all the isospin multi- 
states of (sqsq) exist, respectively, and comparing with 
the regular meson structure of /o (980) (ss). This is help- 
ful in the exotic hadron study in experiments like RHIC. 
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